Consistent model of magnetism in ferropnictides 
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Abstract 

The discovery of superconductivity in LaFeAsO introduced the ferropnictides as a major new class of superconducting 
compounds with critical temperatures second only to cuprates. The presence of magnetic iron makes ferropnictides radically 
different from cuprates. Antiferromagnetism of the parent compounds strongly suggests that superconductivity and magnetism 
are closely related. However, the character of magnetic interactions and spin fluctuations in ferropnictides, in spite of vigorous 
efforts, has until now resisted understanding within any conventional model of magnetism. Here we show that the most 
puzzling features can be naturally reconciled within a rather simple effective spin model with biquadratic interactions, which 
is consistent with electronic structure calculations. By going beyond the Heisenberg model, this description explains numerous 
experimentally observed properties, including the peculiarities of the spin wave spectrum, thin domain walls, crossover from 
first to second order phase transition under doping in some compounds, and offers new insight in the occurrence of the nematic 
phase above the antiferromagnetic phase transition. 



I. INTRODUCTION 



The recent discovery of superconductivity in ferropnic- 
tides has ended the hegemony of cuprates and ushered in 
the new "iron era" in the studies of high-temperature 
superconductivity [H-Q. All parent ferropnictide com- 
pounds, unlike cuprates, are metallic and contain iron 
as the indispensable element. The occurrence of antifer- 
romagnetism in the parent compounds of both cuprates 
and ferropnictides suggests that magnetism and super- 
conductivity may be interrelated, and that the pairing 
mechanisms involve spin fluctuations in both cuprates 
and ferropnictides. However, magnetism in ferropnic- 
tides demonstrates a number of unique and puzzling 
properties, which are inconsistent with any conventional 
model. Understanding the character of magnetic inter- 
actions and spin fluctuations in ferropnictides is there- 
fore of paramount importance for the problem of high- 
temperature superconductivity. In particular, a reason- 
able effective spin Hamiltonian supported by experimen- 
tal measurements and theoretical calculations is highly 
desirable. Until now, such spin Hamiltonian has not been 
found. The purpose of this article is to solve this prob- 
lem. 

The antiferromagnetic (AFM) stripe ground state may 
be obtained in the Heisenberg model with appropriately 
chosen exchange parameters, such as the J\-J-i model 
with J2 > Ji/2, which has been extensively studied in the 
context of frustrated magnetism However, this model 
neglects the large anisotropy of the nearest-neighbor ex- 
change constant, which is firmly established both from 
first-principles calculations [5[ and by neutron scattering 
measurements 0, El 0] ■ The three-dimensional character 
of magnetic interactions is also important @. Without 
these features, the Heisenberg model can not reproduce 
the qualitative features of the spin-wave spectrum. It 
is common to introduce the anisotropy of the exchange 
parameters explicitly; in particular, J\ is split into J\ a 



and J\b for the pairs with antiparallel and parallel spins, 
respectively. Theoretically these parameters can be cal- 
culated using the linear response (LR) technique. The 
anisotropy of J\ strongly depends on the local moment /x; 
Jib (between parallel neighboring spins) can even change 
sign at large values of [i 0, [l(J . Such strong dependence 
of the exchange parameters on the local moment clearly 
rules out the traditional Heisenberg model description for 
ferropnictides. Moreover, the anisotropy of exchange pa- 
rameters breaks the symmetry of the spin Hamiltonian, 
which thereby becomes unsuitable beyond the spin-wave 
regime. In particular, magnetic phase transitions can not 
be described with such anisotropic Hamiltonian. 

Another problem of the Heisenberg model is that it 
fails to describe the energies of noncollinear structures 
connecting the degenerate AFM domains. Specifically, 
in this model with interaction of any range the AFM 
domains with ordering vectors (tt, 0) and (0, it) are con- 
nected by a continuous degenerate set of noncollinear 
states (at least at the mean- field level). In clear con- 
tradiction with this model, band structure calculations 
reveal the presence of a rather high energy barrier [i| Hll | 
between (it, 0) and (0, it) states. In addition, it is difficult 
to reconcile the Heisenberg model with the fact that the 
domain walls between different AFM domains (both twin 
and antiphase) are very thin }12| . 

Failure of the Heisenberg model in ferropnictides led to 
the belief that the spin Hamiltonian approach has to be 
abandoned altogether in favor of a model that explicitly 
takes into account fermionic degrees of freedom. Strong 
anisotropy of the magnetic interaction in the stripe phase 
was argued to be linked with orbital ordering [lj§]. In 
the model with a non-interacting fermion sector added 
to the local moment subsystem, the double exchange is 
anisotropic in the stripe phase [14[ . Although this model 
may be able to reproduce the spin wave spectrum, it con- 
tradicts transport measurements showing a rather mod- 
erate in-plane anisotropy of the resistivity, whose sign is 
opposite to that expected in the double exchange picture 
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15]. In addition, it is not well suited for the studies of 
magnetic thermodynamics. 

Thus, the Heisenberg model of any range clearly fails to 
describe the magnetism of ferropnictides. It is therefore 
natural to request a more refined theory which would (1) 
preserve the symmetry of the lattice, (2) generate the cor- 
rect spin wave spectrum and anisotropic linear response 
properties in the stripe phase, (3) eliminate the spuri- 
ous continuous ground state degeneracy, (4) be suitable 
for the studies of magnetic thermodynamics at finite tem- 
peratures, and (5) be consistent with electronic structure 
calculations. It is also desirable for the model to be ca- 
pable of reproducing subtle features of magnetic thermo- 
dynamics observed in ferropnictides, including the occur- 
rence of both first- and second- order phase transitions 
and the possibility of the nematic symmetry breaking. 
In this article we show that all these requirements can 
be met by including non- Heisenberg interaction in the ef- 
fective spin Hamiltonian. The resulting effective Hamil- 
tonian is found to explain numerous experimentally ob- 
served magnetic properties. 

II. RESULTS 

In the following, we justify the choice of the spin Hamil- 
tonian, present an alternative interpretation of neutron 
scattering experiments, and then describe the magnetic 
thermodynamics in ferropnictides using the mean-field 
approximation and classical Monte Carlo simulations. 

Model spin Hamiltonian 

Our treatment is based on the following effective clas- 
sical spin Hamiltonian: 

H = J2 [jijSiSj - Kn (S 4 S,) 2 ] (1) 

i<j 

The summation in (1) is taken over distinct pairs of lat- 
tice sites. The first term is the conventional (isotropic) 
exchange interaction, and the second term is the pairwisc 
biquadratic interaction, which is the simplest and most 
natural form of non-Heisenberg coupling. We found that 
a consistent description of magnetism in ferropnictides 
is obtained using the J\-K-J%-J c model, which includes 
the nearest and next-nearest in-plane exchange param- 
eters Ji and J2, the interplane exchange parameter J c , 
and the nearest-neighbor (in-plane) biquadratic coupling 

K = ks 2 . 

Note that a weak effective biquadratic interaction 
appears in the J1-J2 model through the "order-from- 
disorder" mechanism Q, and for pnictides it has been 
included in Ref. (l6j assuming a very small interlayer 
coupling. Here we argue that a fairly strong biquadratic 
interaction is already present in the microscopic spin 
Hamiltonian in ferropnictides. Although magnetostruc- 
tural coupling can also effectively induce non-Heisenberg 



magnetic interactions, direct calculations show that 
this distortion in the AFM phase has almost no effect on 
magnetic interactions. 

The presence of non-Heisenberg terms in the Hamil- 
tonian radically changes the interpretation of the LR 
exchange parameters Jjf*, which are defined as second 
derivatives of the total energy with respect to spin rota- 
tions. Rotating the spins at sites i and j in the opposite 
directions by an angle 0/2, we find 

d 2 H 

Jij = — (SjSj) ~Q(p~ = Jij ~ 2KijS e^ej, (2) 

where ej = Si/S. The contribution of the biquadratic 
term depends on whether the spins in the given pair are 
parallel or antiparallel. For the stripe ground state, we 
find J u = Ji + 2KS 2 and J lb = J x - 2KS 2 . Thus, the 
anisotropy of the nearest-neighbor LR exchange param- 
eter in the stripe phase is captured by the biquadratic 
term in Hamiltonian ([1]). 

Spin wave spectrum 

The spin wave spectrum of CaFe2As2 was measured 
using inelastic neutron scattering in Refs. d, 0] and fit- 
ted to the anisotropic J\-J<2,-J c Heisenberg model. Here 
we reinterpret this fitting in terms of our isotropic J%- 
K-J2-J c model, which produces an identical spin wave 
spectrum. The fitted value of K is strongly affected by 
the inclusion of the zone-edge magnons Q • Table 1 shows 
the reinterpreted parameters of this fitting. 
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TABLE I: Parameters of the model (in meV) fitted to spin 
wave spectrum @|. Here K = (Ji a — Jif,)/4, and J s — 4J2 + 
2J C . The parameter J s and Tn are found assuming S = 0.4. 



Fig. Q] illustrates the effect of K on the spin wave spec- 
trum described by the following dispersion relations: 

E(q) = B 2 , (3) 

where 

Aq = 2S[( Ji - 2K) cos(Trfc) + 4K + J c + 2 J 2 ] (4) 
B q = 25 [( Ji + 2K) cos(nh) + 2J 2 cos(nh) cos(?rfc) 
+ J c cos(7rZ)], 

where q is the reduced wave vector (hkl). Starting from 
the values taken from Ref. (Table 1), the evolution 
of the spectrum as K is decreased to zero is shown. As 
noted in Ref. Q, the absence of a minimum at the zone 
edge along the [010] direction signals a sign inversion of 
the parameter Ju,. Using spin Hamiltonian (1), this fea- 
ture is immediately reinterpreted as a signature of a large 
value of biquadratic coupling K > J±/2. 
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FIG. 1: Effect of the biquadratic coupling on the spin wave 
spectrum. The top curves correspond to the model parame- 
ters taken from Table 1. To obtain the other curves, from top 
to bottom, the ratio K/Ji was reduced to 0.4, 0.2, and 0. 



Mean-field phase diagram 

In the presence of biquadratic interaction the single- 
site problem in MFA contains dipole and quadrupole ef- 
fective fields with self-consistency required for both the 
magnetization rrij = (e^) and quadrupole order param- 
eter rriij = (e^ej) — 5y/3. Thus, apart from an AFM 
(stripe) phase with nonzero rrii, an anisotropic quadrupo- 
lar phase may appear with m t = but my ^ (l8| . 

The MFA phase diagram is shown in the inset of Fig. 
[3J The only relevant phase up to K/J% = 3 is the 
AFM (stripe) phase. The most important feature for our 
purposes is the tricritical point at K t = 5J s /24, where 
J s = 4J2 + 2J C . At K — K t the phase transition to the 
paramagnetic phase changes from second order (at lower 
K) to first order. At K < K t the transition temperature 
does not depend on K in MFA, while at K > K t there is 
a gradual crossover to a linear dependence. The tricriti- 
cal point falls well within he range of realistic parameters 
for ferropnictides, and it is further significantly reduced 
by fluctuations (see below). Thus, magnetic Hamilto- 
nian (1) can yield a first-order transition in ferropnictides 
without introducing any magnetostructural coupling (as 
it was done in Refs. [l9(). This is a rather general feature 
of spin Hamiltonians with biquadratic coupling [2(| ■ 

At K w J s there is a triple point where the quadrupo- 
lar phase appears between the AFM and the paramag- 
netic phases. The quadrupolar phase is always sepa- 
rated from the paramagnetic phase by a first-order tran- 
sition, but the AFM-to-quadrupolar transition changes 
from first to second order at a slightly higher K . The 
triple point occurs beyond the reasonable range of K for 
ferropnictides. 



Monte Carlo simulations 

The main results are presented in Fig. [2] We con- 
sidered a few values of the J2/J1 ratio from 0.6 to 
1.25 at fixed J C /J\~Q.2, and a few values of J c /J\ at 
J2/J1 = 0.6. Note that J%j J\ = 0.5 is the stability limit 



of the stripe phase at T = 0, and one can expect par- 
ticularly strong fluctuations in the vicinity of this point. 
The fluctuations are also enhanced at small J c /Ji, where 
two-dimensional behavior sets in. 

As expected, the transition temperature is significantly 
reduced compared to MFA, and this reduction increases 
as J2/J1 approaches 0.5 and as J c /J\ decreases. Inter- 
estingly, the ratio K t /J s demonstrates a similar trend. 
The critical point for the parameters corresponding to 
CaFe2As2 is estimated at 90 K, which, considering the 
simplicity of the model, is in reasonable agreement with 
the experimental range of 140-170 K for 122 compounds. 
Further, this transition is of first order for these param- 
eters, in agreement with experiment. 
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FIG. 2: Transition temperature and order of phase transi- 
tion as a function of biquadratic interaction K obtained by 
Monte Carlo simulations. Temperature is measured in units 
of T^ IFA = J a /3, which is the second-order T N in MFA. Each 
line is labeled by a set of two parameters (J2/J1; J c /Ji). The 
region of first-order transitions (empty symbols) is schemat- 
ically highlighted by shading and thicker lines. The dashed 
line shows the point of inversion of Ju- The inset shows the 
mean-field phase diagram in a larger area of the same param- 
eters. Solid lines 1 and 4 (dashed lines 2, 3 and 5) denote 
second-order (first-order) first transition. A and C are tricrit- 
ical points; B is the triple point. Q denotes the quadrupolar 
phase. 

Experimentally, a second-order AFM phase transition 
in ferropnictides is typically preceded by an orthorhom- 
bic structural transition. Between these transitions the 
so-called nematic phase occurs Q . Theoretically, the ne- 
matic phase is known to occur at J c = Q and was 
predicted to occur at J c <C J\ as well 16]. It is not 
clear whether this phase can survive in the presence of 
substantial interplane and biquadratic interactions. It is 
rather difficult to establish the separation of the nematic 
and AFM phase transitions in Monte Carlo simulations, 
because the AFM correlation length becomes very large 
whenever this separation becomes large, and therefore a 
spurious "quasi-long-range-order" is generated by finite- 
size effects 21]. However, in a certain parameter range 
we have observed a clear signature of the nematic phase 
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in the finite-size behavior of the fourth-order cumulants 
of the nematic and AFM order parameters (their defini- 
tions are given in the Methods section) . Fig. |3K shows the 
separation between the nematic and AFM phase transi- 
tions determined from the crossing points of the fourth- 
order cumulants for these order parameters calculated 
for different D (an example is shown in Fig. [3Jd). The 
separation determined in this way appears to increase at 
larger J2/J1, smaller K, and smaller J c . It never occurs 
in the range of parameters where the phase transition is 
of the first order. 
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FIG. 3: (a) Relative separation of the nematic (T„) and AFM 
(Tjv) transition temperatures as a function of the biquadratic 
interaction parameter K for different values of ]2 = J2/J1 
and j c = Jc/Ji- Open (filled) symbols correspond to the 
first-order (second-order) phase transition, (b) Example of 
the separation of T„ and Tn in the finite-size scaling behav- 
ior of the respective fourth-order cumulants U4 and U4 ; the 
parameters are for the point shown in red in panel a. 



III. DISCUSSION 

Since the K/J ratios are proportional to S 2 , the AFM 
phase transition can change from first to second order 
if S is reduced. Such behavior was indeed observed un- 
der electron doping which reduces the local moments. 
Indeed, while the phase transition in updoped 122 com- 
pounds (fi ~ 0.8 — 0.9/jfj) appears to be first-order (most 
clearly in CaFe9As? \22i )- it changes to second order un- 
der electron doping [23, |24( . While the decrease of ampli- 
tude of the local moments in doped systems has not been 
measured, the critical temperature of magnetic phase 



transition and the order parameter are clearly decreased 
with doping. The band structure studies further support 
this trend. The experimental situation with the local 
moment and order of the phase transition in LaFeAsO 
and other 1111 compounds is not yet fully resolved, but 
recent measurements suggest that they may be closer to 
122 compounds than previously thought |25| . 

In view of the close proximity of the nematic and AFM 
phase transitions and the corresponding limitations of the 
finite-size scaling procedure, the signature of the separa- 
tion of these transitions, while statistically robust, can 
not be viewed as the final proof. Still, it is remarkable 
that the separation appears to survive in an essentially 
three-dimensional case. In addition, contrary to Ref. Ilq . 
the biquadratic interaction is present on the level of mi- 
croscopic Hamiltonian and is not small. Curiously, Fig. 
3 suggests that at certain conditions the separation can 
go through a maximum as a function of K . The origin 
of this behavior remains to be understood. 

Comparison of the effective spin Hamiltonian fitted to 
experiment (Table 1) with the results of first-principles 
calculations leads to an interesting observation. The ex- 
perimental exchange parameters in Table 1 agree quite 
well with LR calculations || [l(| , except that the calcu- 
lated value of K (which can be deduced either from the 
anisotropy of J\ or from the height of the barrier between 
the two stripe states) reaches its experimental magnitude 
only for larger values of the local moment of order 1.4- 
1.8 [Lb- First-principles analysis of the dynamic suscepti- 



bility [171 ] also shows good agreement with experimental 
spin wave spectrum if the local moment is of that order. 
On the other hand, scaling of the spin moment up by a 
factor of two would also scale up the predicted Tn (Table 
1) by a factor of two and bring it in almost quantitative 
agreement with experiment. This observation seems to 
support the notion that the (mean-field) local moments 
in ferropnictides are roughly twice their observed values, 
which may be screened by dynamical fluctuations. 

We have shown that a consistent model of magnetism 
in ferropnictides is obtained by introducing biquadratic 
interaction in the model spin Hamiltonian. This model 
satisfies all reasonable requirements, reproduces many 
experimental observations, and is supported by electronic 
structure calculations. As a word of caution, the suc- 
cess of the localized model spin Hamiltonian should not 
be taken to rule out itinerant effects. On the contrary, 
the effective spin model should be viewed as an appro- 
priate mapping of the complicated itinerant interactions, 
which include complicated rearrangements of the Fermi 
surface and density matrices induced by magnetic order- 
ing. This interpretation is supported by the fact that a 
large biquadratic term appears at relatively low values of 
the local moment. We find it appropriate to speculate 
that strong non-Heisenberg coupling in the effective spin 
model, being a consequence of strong coupling between 
the itinerant electrons and spin fluctuations, explicitly 
suggests spin fluctuations as a possible pairing mecha- 
nism in ferropnictides. 
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IV. METHODS 

We have studied tetragonal lattices D x D x D with 
periodic boundary conditions, where D — 14, 16, 18, 
and 20. Looping over the lattice sites, a new random 
spin direction for the given spin is tried and accepted 
or rejected using the Metropolis prescription. The aver- 
ages are accumulated during a sufficiently long measure- 
ment run after an initial equilibration run. The lengths 
of these runs were adjusted in order to obtain sufficient 
accuracy. The AFM order parameter L is defined as 
L 2 = L\ + L 2 , where L x and L y are the staggered or- 
der parameters with the x and y staggering direction, 
respectively. The nematic order parameter is defined as 
N = | ^ X)i Srm • where the sum is over the 
four in-plane nearest neighbors, and rjij is equal to ±1 
with opposite signs taken for the two perpendicular direc- 
tions. The order of the phase transition was determined 
by analyzing the behavior of the fourth-order energy cu- 
mulant [26| • If this cumulant converges toward the value 
2/3 for all temperatures with increasing D, the transition 
is of second order. On the other hand, if the energy cu- 



mulant develops a minimum at some temperature, which 
becomes sharper with increasing D, the transition is of 
first order 26] . If the transition was found to be of sec- 
ond order, the critical temperatures for the AFM and 
nematic transitions were determined from the finite-size 
scaling behavior of the fourth-order cumulant of the re- 
spective order parameter [271 ] - For first order transitions, 
the transition temperatures for AFM and nematic phase 
transition were found from the peaks of the correspond- 
ing susceptibilities. 
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